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ABSTRACT 

We report detections and constraints for the near infrared Ks band secondary eclipses 
of seven hot-Jupiters using the IRIS2 infrared camera on the Anglo-Australian Tele¬ 
scope. Eclipses in the Ks band for WASP-18b and WASP-36b have been measured 
for the first time. We also present new measurements for the eclipses of WASP-4b, 
WASP-5b, and WASP-46b, as well as upper limits for the eclipse depths of WASP-2b 
and WASP-76b. In particular, two full eclipses of WASP-46b were observed, allow¬ 
ing us to demonstrate the repeatability of our observations via independent analyses 
on each eclipse. Significant numbers of eclipse depths for hot-Jupiters have now been 
measured in both Ks and the four Spitzer IRAC bandpasses. We discuss these mea¬ 
surements in the context of the broadband colours and brightness temperatures of the 
hot-Jupiter atmosphere distribution. Specifically, we re-examine the proposed temper¬ 
ature dichotomy between the most irradiated, and mildly irradiated planets. We find 
no evidence for multiple clusters in the brightness temperature equilibrium tempera¬ 
ture distributions in any of these bandpasses, suggesting a continuous distribution of 
heat re-emission and circulation characteristics for these planets. 
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1 INTRODUCTION 


A secondary eclipse occurs when the emergent flux from 
a planet is blocked by its host star, which allows the di¬ 
rect measurement of the infrared day side temperature of 
a hot-Jupiter. Assembling a sample of secondary eclipses 
is one way to comparatively study the atmospheres of the 


hot-Jupiter population. For example, Cowan & Agol (2011) 


have compared measured effective temperatures to the in¬ 
cident stellar irradiation received by hot-Jupiters, and pro¬ 
posed that the most irradiated planets have heat re-emission 
properties that are different from those for cooler plan¬ 


ets. Updated results from Schwartz & Cowan (2015) have 


constructed ensemble emission spectra of the hot-Jupiter 
population, and placed further constraints on the albedo 


* Based on observations obtained at the Anglo-Australian Tele¬ 
scope, Siding Spring, Australia, 
t E-mail: george .zhouOanu. edu. au 


and heat recirculation efficiences of the most well charac¬ 


terised hot-Jupiters. Triaud (2014) and Triaud et al. (2014) 


have used multi-band eclipse measurements to create colour- 
magnitudes diagrams of hot-Jupiters, enabling an empiri¬ 
cal comparison with the broadband spectral features of M- 
dwarfs and brown dwarfs. The presence of proposed atmo¬ 
spheric thermal inversion features have also been linked to 
stellar activity ( Knutson et al.||2010 ) via comparative stud- 


The majority of secondary eclipse measurements have 
been made using the Spitzer Space Telescope in the IRAC 


bands (e.g. Deming et al. 2014 Shporer et al. 2014 Lan- 


otte et al.|2014). Forty eight planets have now been sam¬ 


pled in eclipse at the IRAC 4.5 /im band. In comparison, 
only 26 planets have been observed in the Ks band from 
the ground. However, ground-based eclipse observations can 
probe shorter wavelengths that are inaccessible to Spitzer, 
and so probe deeper into the planetary atmosphere and ex¬ 
amine different regimes of atmospheric circulation. They 
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also provide an extended wavelength baseline for compar¬ 
ison with spectral models. Despite the challenges inherent 
in ground-based eclipse measurements for hot-Jupiters, a 
number of facilities are now consistently delivering eclipse 


depth measurements. These include the GFHT (Groll et al. 
2010a|b 2011| Wang et alj2013 Croll et al.|2015 1, 200-inch 

at Palomar (e.g. Zhao et al.[[2012b|a 

O’Rourke et al.[[2014 

Shporer et al. 

2014 Zhao et al. [2014 

), and the ESO 2.2-m 

telescope (e.g. 

Ghen et al.[[2014b|a|c). 

ed the series of eclipse 

In Zhou et al. (2014), we introduc 


observations we are performing at the Anglo-Australian 
Telescope (AAT), with the aim of measuring Ks band 
eclipses for a large number of hot-Jupiters in the south¬ 
ern hemisphere. In this paper, we report eclipse measure¬ 
ments and constraints for seven hot-Jupiters: WASP-2b, - 
4b, -5b, -18b, -36b, -46b, and -76b. The observations and 
analysis are described in Section results and comparisons 
between previous observations reported in Section Sec¬ 
tion discusses the eclipse observations in the context of 
the hot-Jupiter colour-magnitude, colour-colour diagrams, 
and brightness temperature - equilibrium temperature rela¬ 
tionships. 


2 OBSERVATIONS AND ANALYSIS 


2.1 Observing strategy and data reduction 

These eclipse observations were performed using the IRIS2 
instrument ( Tinney et al.|[2QQ4 ) on the 3.9-m A AT at Sid¬ 
ing Spring Observatory, Australia. IRIS2 is a IK x IK in¬ 
frared camera with a HAWAII-1 HgCdTe infrared detector, 
read out over four quadrants in double-read mode. The in¬ 
strument has a field of view of 7b 7 x 7b 7 and plate scale of 
0.4486”/pixel. The observing strategy for each eclipse obser¬ 
vation is similar to that described in Zhou et al. (2014): the 


telescope is defocused to broaden the stellar point-spread 
function, reducing the effect of intra- and inter-pixel sys- 
tematics, and preventing saturation of the target and key 
reference stars. Exposure times are set such that the target 
and key reference stars are kept below peak counts of 20,000 
ADU, so as to keep within the regime where detector non¬ 
linearity is minimised (non-linearity of > 1% occurs above 
40,000 ADU). We apply a non-linearity correction to each 
image, following the IRIS2 manual’s prescription, to correct 
the very small (< 0.05%) non-linearity present below 20,000 
ADU. Typical exposure times are < 10 s per exposure. The 
WASP-18 observations were taken using 5s exposures, with 
20 exposures being averaged and saved as a single frame. 
This observing mode was tested to reduced the data vol¬ 
ume of the observations. However we found that it led to a 
reduction in the precision of eclipse timing and depth mea¬ 
surements, and it was not used for any subsequent observa¬ 
tions. For all other observations, single exposures are saved 
and used in the analysis. Observations of six offset positions 
are taken before and after each eclipse sequence to sam¬ 
ple the sky background. These observations are used as flat 
fields in the reduction process. The field is carefully cen¬ 
tred such that the target and key reference stars do not fall 
on bad pixels. The eclipse sequence is performed in stare 
mode, with the telescope guided to minimise drift of the 
field. Dark frames of the same exposure time as the eclipse 


observations are taken before and after each night. Tests of 
darks taken through an experimental night showed no drifts 
in the dark current. Information on the systems observed, 
and the specifics of each eclipse observation, including date, 
number of exposures, median cadence, median point-spread 
function full width at half maximum (FWHM), are given in 
Tabled] 


Each object frame is dark subtracted and flat divided. 
A master dark frame is median combined from darks of the 
same exposure time, taken on the same night as the object 
frames. Master flat fields are created from the set of off¬ 
set frames taken before and after the eclipse sequence, with 
stars masked, and median combined. Bad pixels in the ob¬ 
ject frames are then interpolated over using the surrounding 
pixels via a radial basis function interpolation. Baryocentric 
Julian Date (BJD) time stamps for each frame are calcu¬ 


lated using the converttime task in VARTOOLS ( [Hartman 
et al.||2008 Eastman et al.||2010 ) 


For each frame, stars are identified using Source Extrac¬ 
tor (Bertin & Arnouts 1996) and cross matched using gr- 
match task in FITSH ( Pal|2012 ). Coordinates for the target 
and reference stars are transformed using grtrans, and aper¬ 
ture photometry performed using fiphot For each set of ob¬ 
servations, we extract the photometry through a set of fixed 
apertures. Background flux beneath each aperture is mea¬ 
sured in a 5 pixel wide annulus around each central aperture. 
Adjacent stars within this annular background aperture are 
masked before the background is calculated. We also tried 
extractions using variable apertures, by setting the aper¬ 
ture size per frame as a multiple of the average FWHM of 
the stellar point-spread functions. However, extracting pho¬ 
tometry using fixed apertures yielded light curves with the 
least out-of-eclipse scatter. We conclude photometry using 
variable aperture sizes introduces noise associated with the 
FWHM estimate, and especially for fields with just a few 
bright reference stars. 


2.2 Light curve analysis and modelling 


The object light curve is corrected by a master refer¬ 
ence light curve, constructed from selected reference stars. 
Weights are applied to each reference star light curve such 
that the out-of-eclipse scatter of the object light curve is 
minimised after correction. The final light curve is found by 
minimising the out-of-transit scatter of the object light curve 
for all permutations of object and reference star extraction 
apertures. 

The eclipse light curves are htted with the [Nelson <S| 


Davis (1972) model, using an adapted implementation of 
the JKTEBOP code (Popper & Etzel 1981| Southworth 


et al.|20(M ). The eclipses are modelled with free parameters 


ecosct;, which determines the phase of the eclipse, and the 
surface brightness ratio Sp/Si,, which determines the depth 
of the eclipse. To propagate the uncertainties in the system 
parameters, we also incorporate the free parameters period 
P, primary transit reference time To, planet-star radius ra¬ 
tio Rp/Ri,, normalised orbit radius {Rp + P*)/a, and line- 
of-sight inclination z, each constrained tightly by Gaussian 
priors adapted from the literature uncertainties. In the cases 
where no clear eclipse is seen, we constrain the ecoscj with 
Gaussian priors using literature eccentricity values from pre¬ 
vious Spitzer eclipse measurements, or radial velocity con- 
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Table 1. Targets and observation details 


Target 

ATmag® 

Observation 

Date & Time (UT)^ 

Number of 

Exposures 

Cadence (s) 

Median 

FWHM (pix) 

No. Ref 

Stars 

WASP-2 

9.6 

2014-09-10 08:54-14:46 

3100 

6 

7.6 

9 

WASP-4 

10.7 

2014-09-04 15:54-17:02 

354 

11 

17.4 

4 



2014-09-11 09:23-11:38 

689 

11 

14.2 

4 

WASP-5 

10.6 

2014-09-14 13:17-19:17 

1727 

11 

12.3 

6 

WASP-18 

8.1 

2014-09-05 11:16-14:54 

95 

120^ 

19.0 

3 

WASP-36 

11.3 

2015-03-09 10:21-15:51 

4358 

4 

4.1 

9 

WASP-46 

11.4 

2014-09-11 12:07-17:28 

1600 

11 

13.4 

6 



2014-09-14 08:55-12:53 

1230 

11 

11.6 

6 

WASP-76 

8.2 

2014-09-13 13:26-19:16 

2953 

6 

17.4 

2 


“ 2MASS magnitudes ( [Skrutskie et al.|2006| >. 

^Start and end of each observation sequence. 

^20 X 5s exposures are averaged and used for analysis. Individual exposures were not saved. 


straints. The ecosu constraints, when applied, are noted in 
Table The light travel time has been accounted for when 
htting for the eclipse timing. 

To incorporate the inhuence of instrumental and atmo¬ 
spheric variations into the reported uncertainties, we model 
the light curve as a function of external parameters simul¬ 
taneously with the model ht. This is even more impor¬ 
tant for infrared light curves than at optical wavelengths, 
given the greater pixel-to-pixel sensitivity variations in in¬ 
frared detectors, combined with the much larger variability 
of the infrared sky background. The inhuence of the exter¬ 
nal parameters is modelled as a linear combination of fac¬ 
tors which include (as per Zhou et ^[2014 ): time t, target 
star pixel positions X and y, stellar point-spread functions 
FWHM y, background counts B, and airmass A. Fits are 
performed using all combinations of the external parame¬ 
ters. We then adopt the model with the set of external pa¬ 
rameters that minimise the Bayesean Information Criterion 
(BIG) post-htting. The external parameter model compo¬ 
nents adopted for each eclipse observation are listed in Ta¬ 
ble E] We note in Section |3] when the next-best decorrelation 
models, ranked by BIC, exhibit a different eclipse ht. In the 
cases where eclipse observations are combined from multi¬ 
ple nights (WASP-4b and WASP-46b), the external param¬ 
eters and coefhcients are independent for each night. This 
allows us to account for the different factors that affect the 
observations each time. For example, the seeing conditions 
were stable for the WASP-46b observation on 2014-09-11, 
and variable for 2014-09-14, the instrument model thus con¬ 
tained only the time and airmass components for the hrst 
night, and time, airmass, and FWHM components for the 
second. 

The best ht model parameters and associated un¬ 
certainties are derived using an Markov chain Monte 
Carlo (MCMC) analysis, via the emcee ensemble sampler 
( Foreman-Mackey et al.|2013 ). The MCMC analyses are run 
twice, with the walkers for the second run originating from 
the best ht parameters from the hrst run. In the second run, 
the per-point photon errors are inhated to force a reduced 
= 1. This allows error sources other than photon-noise 


to be included in the uncertainty estimate, and is particu¬ 
larly important for light curves with substantial red noise 
components. 


3 RESULTS 

We detect the eclipses of WASP-4b, -5b, -18b, -36b, and - 
46b at > 3cr signihcance, and provide the 3(7 upper limits for 
the eclipses of WASP-2b and -76b. The full set of derived 
parameters, including hux ratios, eccentricity constraints, 
and brightness temperatures, are listed in TableThe light 
curves for the eclipse observations are plotted in Figures[TJj^ 
with the 3cr upper limits marked where appropriate. 

For WASP-2b, we determine a 3cr upper limit to the 
eclipse depth of < 0.07% from a single eclipse observation. 
The eclipse of WASP-2b has previously been measured with 
Spitzer ( Wheatley et al.|[2010 ) at 3.6, 4.5, 5.9, and 8.0//m. 
Given the lower hux ratio expected in the Ks band com¬ 
pared to the Spitzer bands, and the shallow 0.083 ± 0.035% 
eclipse at 3.6 /xm, our upper limit is consistent with the 
previous observations. The eclipse phase in our model ht 
was constrained by a Gaussian prior centred on the Spitzer 
eclipse detection. The AAT-IRIS2 light curve for the eclipse 
event of WASP-2b is plotted in Figure 

For WASP-4b, we use two partial eclipses to determine 
an eclipse depth of 0.161 q q 4 %, and phase consistent with 
circular orbit of ecosct; = —O.OOl^Q QQg. The Ks band eclipse 
has been previously measured by Caceres et al. (2011), using 
ISAAC on the VLT, at a depth oF 0T85T^^^^%, consistent 
with our measurement to within la. The eclipses were also 


measured by Spitzer at 3.6/xm and 4.5/xm (Beerer et al. 


2011 ). The eclipse phase we measure is also consistent to la 
with that measured by [Caceres et al.| (2011 ) and Beerer et al. 
(2011). The AAT-IRIS2 light curves for the eclipse events of 


WASP-4b are plotted in Figure 

For WASP-5b, we observe an eclipse with a depth of 
0.201 q o 2%5 and an eccentricity estimate from the eclipse 
phase of ecosu = 0-0081q qq 2 . A Ks band eclipse de pth of 
0.269lQ Qg 2 % was measured by 


Chen et al, 


(2014a) using 


GROND on the MPG 2.2m telescope. The eclipse has also 
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Figure 1. Top: The Ks band relative photometry light curve for the eclipse event of WASP-2b as measured on 2014-09-10. Photometry 
from individual observations are plotted in gray, 10 minute bins in black. For each bin containing n points, the error bars are plotted 
to represent the mean per-point uncertainties, which are photon errors inflated to force a reduced = I? scaled by Xj^n. The best 
fit model (eclipse and instrumental parameters) is plotted in red, and 68% of the allowed models reside within the shaded pink region. 
Middle: As above, but with only the 10 minute binned data and the instrument model subtracted. The shaded regions show the allowed 
models in terms of the transit parameters only. Since we can only place an upper limit on the eclipse of WASP-2b, the dashed brown 
line represents the 3(7 regime of the models. Bottom: As above, but showing the data residuals to the best fit model. The shaded region 
shows the difference between each allowed model (eclipse and instrumental) to the best fit model. 


been measured in the J band (Chen et al. 2014a), and at 
the 3.6 and 4.5 fim Spitzer bands ( Baskin et al.||2013 ). Our 
ecoscj measurement indicates the orbit is eccentric with a 
statistical signficiance of 4a, and is consistent within la with 


Chen et al. (2014a), and within 2a with Baskin et al. (2013). 


WASP-5b has a short period, leading to a short tidal circu- 
larisatin timescale. Following Equation 1 of |Dobbs-Dixon| 
et al. ( 20Q4| , the tidal circularisation timescale should be 


1 Myr (assuming a tidal quality factor of Q' = 10^). It is 
interesting that the eccentricty is non-zero, suggesting that 
perhaps a larger tidal quality factor is required to describe 


the system. The AAT-IRIS2 light curve for the eclipse event 
of WASP-5b is plotted in Figure 


For WASP-18b, we measure an eclipse depth of 
0.141 q q 3 %, and an eccentricity estimate from the eclipse 
phase of ecoscj = 0.012lQ QQg. The eclipses have previously 
been measured by Spitzer at 3.6 and 4.5/xm by |Nymeyer| 


et al. (2011); Maxted et al. (2013), and at 5.9 and 8.0/xm 
by Nymeyer et al. I ( |2011| |. The eclipse phase we derive is 
consistent (at the 2a level) with that expected from circular 
orbit, and the measurements from the Spitzer observations. 
The long cadence observing strategy used for the WASP-18b 
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Figure 2. Light curves for the eclipse events of WASP-4b as observed on 2014-09-04 and 2014-09-11, labelled as per Figure[^ but with 
the light curves from each night arbitrary offset for clarity in the top panel. 
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Phase 

Figure 3. Light curve for the eclipse event of WASP-5b, plotted as per the description for Figure 


observations results in a larger uncertainty in eclipse phase 
than for other targets. We also caution that the second-best 
model, ranked by A BIG (with A BIG of 4), involves the 
time and background flux terms, reducing the surface flux 
ratio to Sp/S^, = O.lOlo'oe- The third-best detrending model 
involves the detector y-position, and gave an eclipse depth 
consistent with that of the best fit model. The AAT-IRIS2 
light curve for the eclipse event of WASP-18b is plotted in 
Figure 

For WASP-36b, we report the first eclipse detection 
of the system at depth of 0.131 q q 4 %, and eccentricity of 
0.0041 q qq 5 , consistent with a circular orbit. No previous 
eclipse observations have been reported for this planet. The 
AAT-IRIS2 light curve for the eclipse event of WASP-36b 
are plotted in Figure 

For WASP-46b, we measure an eclipse depth of 
0.261q q 3 %, and ecosu = 0.0041q;qq4, as measured from two 


full eclipse observations. The Ks band eclipse has previously 
been measured by Ghen et al. (2014c), with a reported depth 
of 0.253l°;°“9 


consistent with our measurement to better 


than 1(7. Ghen et al. (2014c) also measured the J and H 


band eclipses using GROND on the ESO 2.2m. The eclipse 
timing is also consistent with that of a circular orbit, and 


that reported by Ghen et al. (2014c) to Icr. The AAT-IRIS2 
light curves for the eclipse events of WASP-46b are plotted 
in Figure 

For WASP-76b, we find a 3cr upper limit eclipse depth 
of 0.3%, with a marginal detection at 2.3a of a 0.13lo‘o6^ 
eclipse. The eclipse phase fit was constrained by a Gaus¬ 
sian prior on the orbit eccentricity ( West et al.||2013" ). The 
marginal eclipse detection is not well constrained in eclipse 
phase due to the lack of sufficient pre-ingress baseline. As a 
result, the 3a upper limit is asymmetric about the eclipse 
centre. No previous eclipse observations have been reported 


























































Secondary eclipses from the A AT 7 



for this planet. The AAT-IRIS2 light curve for the eclipse 
event of WASP-76b is plotted in Figure 

For each set of eclipses, we also calculate a 13 factor 
to check for any residual time-correlated noise to the light 
curves (Winn et ah 2008). We compared the progressive 
binned scatter of the light curve residual (i.e. the data with 
the eclipse and external parameter models subtracted) with 
the expected scatter assuming only photon noise. For every 
m bins of n points, we measure a root mean square (rms) 
scatter (Jn, which is then compared to the expected photon 
noise scaled rms from the unbinned light curve (cri), accord¬ 
ing to: 


tween 60s and 600s, along with the rms-bin size relation¬ 
ships, are shown in Figure 


C7n 



( 1 ) 


For light curves with no time correlated noise, /3 = 1. The 
average (3 value for each light curve, calculated for bins be- 
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Figure 5. Light curve for the eclipse event of WASP-36b, plotted as per the description for Figure 
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Figure 6. Light curves for the eclipse events of WASP-46b, as measured from 2014-09-11 and 2014-09-14, plotted as per the description 
for Figure 
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Figure 7. Light curve for the eclipse event of WASP-76b, plotted as per the description for Figure [^ 
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Figure 8. The rms of the light curve residuals, with the eclipse and external parameter models subtracted, as a function of the bin sizes. 
The red line plots the expected rms assuming no time correlated noise, where the binned scatter scales 1/v^- We calculate the average 
(3 factor for each set of observations, which measures level of time correlation within the data binned at timescales between 60s and 600s. 
/3 > 1 indicates the presence of time-correlated noise in the residuals, while /3 = 1 indicates residuals scale with the bin size. 












Table 2. Model fit and derived eclipse parameters 
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4 DISCUSSION 

We report new Ks band eclipse depth measurements and 
constraints for seven hot-Jupiters. These results bring the 
total number of planets with eclipses monitored at 2 . 1 /xm 
to 25, a sample large enough to allow some initial statistical 
insight into the atmospheres of the hot-Jupiter population. 

Such statistical analyses require a set of robust measure¬ 
ments with reliable uncertainty estimates. To characterise 
the robustness and repeatabilitiy of our measurements, we 
performed an independent analysis for the two full eclipses 
of WASP-46b that we obtained (Figure [^. We derived 
self-consistent eclipse phases and depths between the two 
eclipses to within la (ecoscj of 0 . 0061 q qq 5 and O.OOIIq qq^, 
and eclipse depth F^jF^ of 0.251 q q 5 % and 0.321 q;q 5%). On 
both nights, the eclipses were detected at > 5cr significance, 
and the measured eclipse parameters were consistent with 
each other at the la level. Literature Ks band eclipse mea¬ 
surements also exist for WASP-4b, -5b, and -46b. The eclipse 
depths we report are also consistent within la to all the pre¬ 
vious measurements. This increases our confidence of the un- 
certaintiy measurements presented by the series of ground- 
based eclipse observations to date. For the wider sample 
of literature eclipse measurements, the uncertainties in the 
eclipse depths are also often underestimated: the scatter in 
repeated eclipse depths reported for the same planet is 1.4 
times larger than the mean error estimates for Xs-band ob¬ 


servations (Zhou et al. 2014), and two times larger for Spitzer 
measurements ( Hansen et al.||2014 ). 

However, eclipse measurements remain intrinsically dif¬ 
ficult due to the low signal-to-noise nature of the planetary 
eclipse and the variety of systematic signals that can be in¬ 
troduced. Time-correlated noise often still remains in our 
observations despite decorrelation against the instrumental 
model, as most of the light curve residuals have /3 > 1. In 
cases where the time-correlated noise is significant (e.g. ob¬ 
servations of WASP-36b and WASP-76b), we notice the un¬ 
certainties in the eclipse depth are larger than other eclipses 
of equivalent depths. 

With these cavaets, we can place our observations in the 
context of the hot-Jupiter atmospheres sample. In this sec¬ 
tion, we empirically examine the colour-magnitude distribu¬ 
tions and brightness-equilibrium temperature distributions 
of the hot-Jupiter population. 


4.1 Colour-magnitude diagram 

We compare the broadband colours and magnitudes of the 
planets examined in this study with other hot-Jupiters, 
brown dwarfs, and late M-dwarfs. Brown dwarfs and late M- 
dwarfs have similar effective temperatures as the equilibrium 
temperatures of hot-Jupiters, but higher surface gravities, 
and experience different levels of irradiation to hot-Jupiters. 
In addition, brown dwarf atmospheres are relatively better 
understood, and brown dwarf spectral models are the source 
from which most hot-Jupiter atmosphere models are built. 
Colour-magnitude and colour-colour diagrams can help com¬ 
pare the broadband spectra hot-Jupiters, and examine for 
differences between the atmospheres of hot-Jupiters and 
brown dwarfs (e.g. Marois et al.||2008 Triaud||2014 Triaud 
et al.||2014| ). 

To calculate the absolute magnitudes of the hot- 


Jupiters, we first derive absolute magnitudes for their host 
stars. The published effective temperature, surface gravity, 
and metallicity of the host stars are fitted to the Dart¬ 
mouth isochrones ( Potter et al.|[2008 ), from which the ab¬ 
solute magnitudes at iLs and the Spitzer IRAC bands are 
extracted. To propagate the uncertainties, we draw 10^ iter¬ 
ations of the stellar atmospheric parameters from Gaussian 
distributions about their reported mean and uncertainty val¬ 
ues. For the seven host stars that have Hipparcos distances 
(identified in Triaud||2014 ), the isochrone derived J, if, K 
distance modulus agree with that measured from parallax, 
with 0.2 mag scatter in the residuals. We set 0.2 mag as 
the lower limit for the absolute magnitude uncertainties we 
calculate from isochrones. We then calculate the absolute 
magnitudes of the planets using the measured eclipse depth 
values and uncertainties. In the cases where repeated obser¬ 
vations are available, we calculate a weighted mean eclipse 
depth and a weighted standard error in the mean for the 
uncertainty. The literature eclipse measurements are listed 
in Appendix [A a nd gathered part ially from the Ex oplanet 
Orbit Database^ ( Han et al.|20I4 ) and Table 3 from Bailey 
( 20I4| ). Hot-Jupiters in high eccentricity orbits (HD 80606b, 
WASP- 8 b) were excluded from the list as they are not rep¬ 
resentative of the hot-Jupiter sample. 

Figurej^plots the iLs — [3.6] and iFs — [4.5] colours of the 
hot-Jupiters and brown dwarfs against the M[ 3 . 6 ] and M[ 4 . 5 ] 
band absolute magnitudes. These bands are chosen as they 
have the most number of eclipse measurements. The hot- 
Jupiters are plotted in colour to represent their equilibrium 
temperatures, with point sizes indicating their planet radii. 


Brown dwarfs compiled from Dupuy & Liu (2012) are plot¬ 


ted in gray scale to represent their spectral classes. Model 
colours from BT-Settl ( Allard et al.| 20 I 2 ), with abundances 
from Asplund et al. (2009), are plotted for reference. 


The sampled hot-Jupiters reside around the M-L spec¬ 
tral classes. The colours of hot-Jupiters are consistent with 
the colours of brown dwarfs. The Ks — [4.5] colours are 
marginally redder for the hot-Jupiters than brown dwarfs. A 
discrepancy between the two populations based on 4.5/xm- 


related colours was suggested by Triaud et al. (2014), with 
the mechanism being the absence of absorbing spectral fea¬ 
tures at 4.5/xm for hot-Jupiters compared to brown dwarfs. 


The transiting irradiated brown dwarf KELT-Ib (Siverd 
et al.|2012 ) is also plotted. It is the only brown dwarf, receiv¬ 
ing similar irradiation levels as hot-Jupiters, with secondar 
eclipse measurements. The colours of KELT-Ib matches well 
with that of isolated brown dwarfs and brown dwarf atmo¬ 
sphere models. 

Eor a direct comparison of the spectral properties of hot- 
Jupiters and brown dwarfs, we plot their iLs — [3.6] ys Ks — 
[4.5] colour-colour relationship in Eigure |I0[ This removes 
the luminosity dependence on radius and reduces the scatter 
in the distribution. Here, it is interesting to note that the 
current sample of hot-Jupiter colours matches well with that 
of brown dwarfs. 

We expect a greater diversity in the hot-Jupiter spec¬ 
tral properties compared to brown dwarfs. Hot-Jupiters have 
radii that differ by a factor of two between similar mass plan- 


^ This research has made use of the Exoplanet Orbit Database 
and the Exoplanet Data Explorer at exoplanets.org. 
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Figure 9. Colour-magnitude diagrams for hot-Jupiters and brown dwarfs. Hot-Jupiters are plotted in colour according to their equilibrium 
temperature, and with point sizes relative to their planet radii. The planets with Ks band eclipses reported in this paper, and with 
available Spitzer eclipse observations, are labelled. Where appropriate, upper limits in the colour axis are given. The transiting brown 
dwarf KELT-lb is plotted as the black diamond. It is the only irradiated brown dwarf with secondary eclipse measurements available. 
The colour-magnitudes of brown dwarfs with parallaxes and photometry compiled by |Dupuy fc Liu| ( |2Q12| > are plotted for comparison, 
the points are gray-scaled according to spectral class. The BT-Settl model colours, from < 4000 K, are marked by the black lines. The 
corresponding effective temperatures of the models are marked by crosses every 500 K. Models at log^f = 3.5 and 5.0, corresponding to 
the surface gravities of Jupiter and typical brown dwarfs, are marked by the dashed and solid lines respectively. The absolute magnitudes 
for the models are converted assuming 1 Rj objects. References for the parallaxes of the points plotte d are|Andrei et al.|(|2011|>; Artigau 


et al. ( 12010 1; [Benedict et al.|(|1999|>;|Costa et ar](|2005||2006 l;|Dahn et ^(|2002]);|Dupuy Sz Liu|(|2012|); [Gatewood Sz Coban|(|2009 | ;|Geye~ 


Harrington et al.[ ( [1993[>; [Henry et al.[([2Q06 >; Kirkpatrick et al.[ ([2011 >; [Lepine et al. (2009 rT p^arocco et al.[ (|2010[ > ; [Monet 

Reid et al.[ ([200^; [Schilbach et al.[ ([2009 1; Subasavage et aL] ~( 2009|~ Teixeira et al. (20081; jTinney et al.[ ( [l995[ >; [Tinney 


et al. 


19881; 


et al. 


19921 


il996 


; [Tinney et al.j ([2003 Altena et al.[ ([1995[>; 


van Leeuwen 


_ ^ _. ._. ._. 

used are [Biller et al.[(|2010[>; Bonnefoy et al.[ ([201 ![>; [Chiu et al.[([2006[>'; [Close et al.[([2002[>;' Currie et al.[([2011[>;[Cutri et al.[([2003[>; Deacon[ 


Vrba et al.[ ([2004[). References for the 2MASS photometry 


et al.[ ([2012J; [Dupuy et al.[ ([2009a|b[ [2010J; [Dupuy & Liu[ ([2012J; 

Esposito et al.[ ([2013 

1; [Forrest et al.[ ([1988J; 

Galicher et al. 

2011 

1; 

>; 

Golimowski et al.[( 

2004J; [Henry & McCarthy[([1993[);[Hewett et al.[( 

2006[);[Janson et al.[( 

2011[);[Jones et al.[([1996 

); [Kasper et al. 

2007 

King et al.[ ( 

to 

o 

o 

Kirkpatrick et al.[( 

201 IJ; [Knapp et al.[ ([2004J; 

Konopacky et al.[([2010 

|;[Lane et al.[ ([2001 [); [Leggett et al.[( 

1998 

[2000[ 

2001[ [2002a 

b[[2007[);[Lm et al.[([2006 

[2008[ 2010[); [Lodieu et al. 

([2007J; [Lowrance et al. 

(|2000|l;|Marois et al.|(|2008l|2010|l; 

Metchev &[ 


Kirkpatrick et arf([2011[>; [Leggett et al. |(|2007|l;|Liu et al.|(|2010Ti; 


; [Reid fc Cruz| l|2002^ . References for the IRAC photometry used are 


Luhman et al.[ ( [2012[ >; [Patten et al.[ ( [2006| ; [Wright et ah] ( [2010[ >. 


ets, and have atmospheres that are heated from above and 
below at different levels of irradiation, resulting in a wide 
range of possible pressure-temperature profiles. We there¬ 
fore expect the colour distribution of hot-Jupiters to exhibit 
significantly greater scatter than that of brown dwarfs. How¬ 
ever, the scatter in the colour distribution is currently dom¬ 
inated by the measurement uncertainties of the eclipse ob¬ 
servations. Further repeated observations, for robust colour 
measurements, and a greater sample size, may help distin¬ 
guish between the brown dwarf and planet population. 


4.2 Brightness — equilibrium temperature 
distribution 

3D models investigating the circulation of strongly irradi¬ 
ated hot-Jupiters have predicted large day-night tempera¬ 
ture differences, with the most irradiated planets developing 
strong, super-rotating, equatorial jets and large longitudi¬ 


nal temperature gradients (e.g. Showman & Guillot 2002 
Dobbs-Dixon Lin|2008 ). These have been revealed by ob¬ 
servations of infrared phase curves with peaks offset from 
the sub-stellar point (e.g. Knutson et al.|2007 Zellem et al. 


2014 Stevenson et al.||20i4b ). Showman et al. (2015) have 

shown the extent of this day-night temperature difference is 
dependent on the level of irradiation and the rotation rate 
of the planet. In comparison, the circulation of mildy irradi¬ 
ated and/or rapidly rotating hot-Jupiters are expected to be 
dominated by latitudinal variation, and weaker longitudinal 
differences. 

To probe for a boundary between the mildy irradiated, 
thermally well mixed hot-Jupiters, and the strongly irra¬ 


diated hot-Jupiters, Cowan & Agol (2011) and Schwartz 
& Cowan (2015) used the available multi-band eclipse ob¬ 


servations to compare day-side effective temperatures of 
hot-Jupiters with their expected equilibrium temperatures. 
They tentatively identified two populations of hot-Jupiters, 
with the most irradiated planets having lower heat recircu- 
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Figure 10. Ks — [3.6] vs Ks — [4.5] colour-colour diagram for hot-Jupiters and brown dwarfs. Colour-colour diagrams reduce the effect 
of the large scatter in planet radii, therefore luminosity, on the spectral class comparisons. The plot markings and references are the 
same as Figure 



Figure 11. The normalised brightness temperature Tg /Teg - equilibrium temperature Teg distribution from secondary eclipses measured 
at near infrared bands. The hot-Jupiters with Ks band eclipses reported in this paper, and with relevant Spitzer band observations, 
are labelled. Gaussian Mixture Model (GMM) clustering preferentially selects a single component model for the distribution at all the 
bands. The blue and red colours show the clusters if we force a N = 2 component fit to the distributions. 
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Figure 12. The p-value - Teq relation from a moving K-S test along the Ts/Teq - Teq distribution. Minima in the p-values should 
indicate divisions between distinct populations in the sample. We find no significant p-value minima (all > 0.05), indicating a lack of 
significant division in the population at all the bands. 


lation efficiencies, and higher relative effective temperature, 
than mildly irradiated hot-Jupiters. 

We re-examine this proposed dichotomy for the pho¬ 
tometric bands where a significant number of hot-Jupiters 
have been sampled in eclipse. Each infrared band probes 
a diherent layer of the planetary atmosphere, with shorter 
wavelengths probing higher pressure regions, where models 
predict better thermal mixing than at higher altitudes. For 
each planet sampled in each of the Ks and Spitzer IRAC 
bands, we calculate an equilibrium temperature Teg, assum¬ 
ing zero albedo and no heat redistribution. The uncertainty 
on the equilibrium temperature is calculated from 10^ it¬ 
erations of random sampling, such that errors in the semi¬ 
major axis, stellar radius, and stellar ehective temperature 
are propagated. For each band, we calculate a brightness 
temperature Tb from the reported average eclipse depths 
(listed in Appendix 0 . The Teq of each planet and each 
band is plotted against the normalised TsjTeq ratio in Fig¬ 
ure [TT] 

To test for multiple populations within the TsjTeq dis¬ 
tribution, we apply the moving two-sample Kolmogorov- 
Smirnov (K-S) test and the Gaussian Mixture Model 
(GMM) analyses. These techniques have previously been 
used to check for distinct populations in the metallicity - 
planet radius ( Buchhave et al.[|2Q14 Schlaufman||2015 ) and 
metallicity - planet period distributions ( Zhu||20l'5 ). 

In a moving K-S test, we split the population into two 
samples along a series of Teq values, and in each instance 
calculate the K-S test p-values to test the null hypothesis 
that the Tg/Teg distribution for the two samples originate 
from the same population. If a significant minimum is ob¬ 
served in the Teq - p-value relationship, then we can state 
that the distribution is made of two distinct populations. To 
propagate the uncertainties in Teq and Tb, we perform the 
moving K-S test 10^ times, at each iteration drawing each 
point from distributions about its mean and error. The Teq 


p-value relationship from the moving K-S test is plotted in 
Figure We recover the weak division suggested by |Gowan| 
Agol| ~i2011 ) in all the Spitzer bands, finding a tentative 


division at 2170, 2590, 2440, and 2380 K in the 3.6, 4.5, 5.8, 
and 8.0/im bands. There appears to be no significant divi¬ 
sion in the Ks band, with the lack of a clear minimum in 
the Teg - p-value relationship. However, the division is sta¬ 
tistically insignificant in all of the bands, with the minimum 
p-value consistently ^ 0.05. 


Schlaufman (2015) suggested modelling the populations 


as Gaussian mixtures as a more rigorous way of distinguish¬ 
ing between multiple clusters in the population. We employ 
the GMM clustering function in the Python package Scikit- 
learn ( Pedregosa et al.||2012 ). We fit the Teg-Tg/Teg distri¬ 
bution with GMMs consisting of N — 1, 2..., 5 full Gaussian 
components. The model that minimises the BIG is chosen 
as the best fit model. To take into account the per point 
uncertainties, we draw the population from their measure¬ 
ment uncertainties 10^ times, each time performing the BIG 
calculation and model selection. We find that the single com¬ 
ponent model is preferred for the Teg-Tg/Teg distribution 
at every band. The N ^ 2 component models are rejected 
> 90% of the time in the Ks, 4.5, 5.8, and 8.0/xm bands, 
and > 60% of the time for the 3.6 /xm band. We find a lack 
of conclusive evidence that the eclipse sample can be split 
into two populations. 

Given the limits of the current data, these tests show 
that the atmospheric circulation properties of the hot- 
Jupiter population is continuous. We suggest that there is 
likely no sharp divide between the warm Jupiters dominated 
by latitudinal circulation and the hot-Jupiters with longitu¬ 
dinal circulation. 

However, if we were to force a two component fit to 
the GMM, we find a relatively consistent result between the 
groups identified in each band. There is also a consistency in 
the 4.5, 5.8, and 8.0/xm bands between the divisions identi- 
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fied by the moving K-S test and the division between clusters 
found by GMM. Figure [^colour-codes each planet accord¬ 
ing to the classifications from the two-component GMM fit. 

The difference between the mean Ts/Teq between the 
two clusters is increasing with wavelength. The difference in 
the mean is 0.03 ± 0.18 at the Ks band, and 0.24 ± 0.20 
at [8.0]. Longer wavelength probe the upper planetary at¬ 
mosphere, where the day-night temperature gradient is ex¬ 
pected to be highest, whilst better thermal mixing is ex¬ 
pected at higher pressures deeper in the atmosphere. The 
same effect should also lead to a larger difference between 
the two proposed hot-Jupiter populations at longer wave¬ 
lengths, as demonstrated by this trend. We suspect that an 
underlying smooth transition may exist, but must await a 
larger sample size before becoming statistically significant. 
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APPENDIX A: LITERATURE ECLIPSE 
MEASUREMENTS 

Table [AT] presents literature eclipse depth measurements for 
the As, Spitzer IRAG 3.6, 4.5, 5.8, 8.0 fim bands. 
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Table Al: Literature secondary eclipses in the Ks and Spitzer IRAC bands 
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U.^i9_o.oi5 

U.19»1_0 0069 

0 401+0.043 
U.9Ui_Q Q 43 

0 24+0-026 
'^•^^-0.026 

KELT-lb 

U.iDU_o 020 

u.iyo_o 010 

U.ZUU_o 012 



Kepler-5b 


0 104+0-017 

U.iUC)_o 017 

0 107+0-015 
-0.015 



Kepler-6b 


u.uon_Q q27 

(4 -,.-,+0.019 

U.lOl-0 019 



Kepler-12b 


0 1 47+0-02 

U.191-0 02 

0 116+®-®^^ 
U.11D_0 031 



Kepler-13b 

n 1 99+0.051 
U.iZZ_o 051 

U.iDD_o 031 

n 999+0.023 

u.zzz_Q q23 



Kepler-17b 


0-25l°:°° 

n ^1 +0.035 
^•'^^-0.035 



OGLE-TR-113b 

n 1 y+O.OS 
-0.05 





Qatar-lb 

0 1 3fl+0-034 
U.lOD_0 034 





TrES-lb 


0.0831^21 

'9.UOO_o 013 

u.unn_Q q24 

U.iOZ_o 042 

n 994+0.036 

u.^^n_o 036 
n 91 q+0.042 

U.Zi9_o 042 

TrES-2b 

' 9 .OOZ _0 oil 

0 127+®-°^^ 
-0.021 

n oq+0.024 
u.^9-0 024 

u.inn_o 054 

0-3591°:°° 

TrES-3b 

U.zni_o 043 
n 1 qq+0.018 

'9.iOO_o 016 

0 346+0-035 
U.9no_0 035 

n Q79+0.054 
U.91Z_o 054 

0 449+0-097 

u.nnn_o 097 

0-4751°:“° 

TrES-4b 


0 147 + 0.011 

U.191-0 oil 

n 140 + 0.016 

0. 148-0 016 

'9.ZDi_o 059 

n Q 1 q+ 0.044 
U.9i9_o 044 

WASP-lb 


O-lS+nie 

0 21 7+0.017 
' -0.017 

0 274+®‘®^® 

u.^/^-O 058 

0 474+0.046 

’0-^'^-0.046 

WASP-2b 

< 0.07 

0.083ini 

u.ion_o 017 

0 102+0-077 
U.iyz_o 077 

0 235+0-059 

'9.Z90_o 059 
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Table A1 - continued from previous page 


Planet 

Ks depth (%) 

[3.6] depth (%) 

[4.5] depth (%) 

[5.8] depth (%) 

WASP-3b 

0.1811°;“ 

0.193l°;“t 

u.zu»_o 028 

0.282l°;“4 


WASP-4b 

u.ioo_Q 03^3 

0 in+O-04 
'0-l'^-0.04 

n.c)iy_o 031 

0.3431^27 


WASP-5b 

QQ 2 

U.ZU_Q Q2 

0 197+®-°^® 

U.j.a < -0.028 

n 997+0.024 

U.ZC)/_o 024 


WASP-lOb 

n 1 97+0-013 

u.io/ -0 019 




WASP- 12 b 

u.z»»_o q65 

U.zyD_o 014 

0.42ll°;°“ 

0.428l°;“4 

0.6961°;°° 

WASP-14b 


0-2241°;“ 

0 224+0-018 
U.ZZ^_0 018 


WASP-17b 



u.zz»_o 013 


WASP-18b 

f) 1O+0.03 
’^•l'^-0.03 

n O+ 0.02 
’^•'^- 0.02 

0 . 304 IHI 9 

n 07 + 0.03 
-0.03 

0 379+0-015 

n.c)/y_o 015 

n 07 + 0.03 
-0.03 

WASP-lOb 

0.3661^72 

0.287l°;“ 

0.4831^25 

n 572+®-®^ 
U.O/Z_o 03 

0.65l°-+ 

WASP-24b 


U.iOy_o 013 

u.zuz_o 018 


WASP-33b 

-0.04 

u.z^^-o 02 

'-'•^0-0.05 

0 41+0-02 
^•^-'-- 0.02 


WASP-36b 

n 1 O+0.04 
'-'•-^'^-0.04 




WASP-43b 

0 194+®-®^® 
u.±»^-0 029 

U.iy/-o 042 
nisi +0-027 
U.i8i-o 027 

0 347+0-013 

u.oa/ -0 013 

n.ooz_o 015 


WASP-46b 

n 9K9+0-063 

n.zoc)_o 06 

U.ZO-o 04 




WASP-48b 

u.iua_o 027 

0.176l°;“J 

0 214“^®'®^ 
n.zi^_o 02 


WASP-76b 

< 0.3 




WASP-80b 


0.0455l°;“ 

U.Uy4:^_0 0065 


XO-lb 


0.086inS^ 

n 1 99+0.009 

U.izZ -0 009 

0.2611^31 

XO- 2 b 


O.OSliniv 

0.0981^2 

0 1 67+0-036 
U.iD/-o 036 

XO-3b 


0 ini+O-004 
U.iUi-o 004 

0.1431^06 

0 1 34+0-049 
U.id^-O 049 

XO-4b 


n.UOD o 006 

n i3>5+®+’^ 

U.iC )0 o 007 



[8.0] depth (%) 
0.328l°;°|8 


u.oyD_Q 096 

fK -| 0-1 +0.022 

U.i8i_o 022 

n 907+0.039 
u.zo/ _o 039 

0 41 + 0-02 
’0-^J--0.02 

0 79+0-12 
O'-''J-0.12 
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